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ABSTRACT 

Expendable  Bathythermograph  (XBT)  temperature  depth  profile  data 
from  17  locations  in  the  north  Pacific  Ocean  has  been  analyzed  to  extract 
information  concerning  the  thermal  microstructure.    The  microstructure 
information  is  then  examined  to  determine  if  the  distribution  can  be 
associated  with  the  climatic  regions  of  the  north  Pacific  Ocean  proposed 
by  Tully  [  1964  ]. 

A  modified  binomial  filter  used  by  Roden  [  1971  ]  to  analyze  STD  data 
was  applied  to  find  a  mean  temperature  series.    A  temperature  per- 
turbation series  results  when  the  mean  is  subtracted  f  rom  the  original 
XBT  series.    Two  resulting  parameters ,  average  mean  temperature 
perturbation  vertical  thickness  scale  and  average  maximum  amplitude 
of  the  temperature  perturbations ,  can  be  associated  with  the  micro- 
structure. 

The  average  mean  temperature  perturbation  vertical  thickness  scale 
was  found  to  increase  with  latitude  from  20°N  to  55°N ,  and  as  the  thick- 
ness scale  increased,  the  maximum  amplitude  of  the  perturbations  was 
found  to  decrease. 
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I.    INTRODUCTION 

Temperature  and  salinity  distributions  in  the  ocean  determine  the 
propagation  paths  for  acoustical  energy.    The  gross  temperature  structure 
is  used  for  ray  path  calculations  to  predict  transmission  characteristics 
[Cohen  and  Weinberg,  1969].    However,  sound  measurements  at  one  location 
reveal  a  fluctuating  acoustic  pressure  level  [Sheehy,  1950]  attributed  to 
variations  in  the  fine  temperature  and  salinity  structure  that  cause  small 
but  important  changes  in  the  propagation  paths.    At  higher  frequencies 
the  microstrucrure  becomes  increasingly  significant.     A  theory  for  acoustic 
wave  transmission  has  been  developed  for  a  medium  containing  thermal 
microstructure  assuming  that  the  microstructure  can  be  represented  by 
inhomogeneities  of  a  thermal  nature  that  are  spherical  or  spheroidal 
"patches"  randomly  distributed  in  space  [For  example,  Liebermann,  1951; 
Mintzer,  1953a,  1953b,  1954].    The  measured  sound  fluctuations  have  been 
related  to  the  focusing  or  defocusing  of  acoustic  waves  by  these  "patches" 
at  ranges  less  than  ka     (where  k  is  the  acoustic  wave  number  and  a  is  the 
effective  radius  of  the  patch)  and  to  interference  at  ranges  greater  than 
ka2  [Skudrzyk,  1957]. 

Recent  findings  concerning  the  microstructure  do  not  support  this 
"patch"  distribution  assumption.     A  layered  microstructure  with  layers 
of  small  gradient  a  few  to  a  few  tens  of  meters  thick  separated  by  sheets 
of  high  gradient  a  few  to  a  few  tens  of  centimeters  thick     id  with 
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horizontal  dimensions  up  to  thousands  of  meters  have  been  found  in  many 
regions  of  the  oceans  [Cooper,  1967;  Cooper  and  Stommel,  1968;  Neshyba, 
Neal,  and  Denner,  1969;  Osborn,  1969;  Roden,  1971;  Simpson,  1971; 
Stommel  and  Fedorov,  1967;  and  Woods  .  1968].    Several  mechanisms  have 
been  proposed  for  generating  the  layered  microstructure ,  and  it  may  be 
that  a  different  mechanism  is  dominant  for  different  ocean  regions. 
This  layered  microstructure  concept  requires  a  re-examination  of  the 
present  theory  for  the  transmission  of  acoustic  energy  in  a  randomly 
inhomogeneous  medium. 

In  1957,  Whitmarsh,  Skudrzyk,  and  Urick  concluded  a  paper  on  the 
scattering  of  sound  in  the  sea  with  the  following  comment:    ". . . .  the 
greatest  g?p  seems  to  lie  in  the  micrcthermal  statistics  of  the  ocean. 
If  these  were  known  better  than  they  are  now,  apparently  the  amplitude 
statistics  of  the  fluctuations  of  sound  transmitted  through  the  sea  could 
be  predicted  roughly,  but  with  some  assurance.  "       Some  11  years  later, 
during  a  meeting  on  Turbulence  In  The  Ocean  held  at  The  University  of 
British  Columbia,  R.  W.  Stewart  "pointed  out  that  a  case  existed  for 
another  degree  of  exploration  of  the  ocean.    While  in  the  past,  widespread 
surveys  have  sought  to  measure  the  temperature  and  salinity  fields ,  we 
must  now  concentrate  on  a  widespread  search  for  fluctuations  in  the 


Whitmarsh,  D.  C.  ,  Skudrzyk,  E.  ,  and  Urick,  R.  J.  ,  "Forward  Scattering 
of  Sound  in  the  Sea  and  Its  Correlation  with  Temperature  Microstruc- 
ture ,"  J_1_Ac^u£._Soc_Am.  ,  v.   29,  p.  1141,  October  1957. 
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2 

temperature  and  salinity  fields.         The  only  data  that  exists  with  such 

broad  coverage  is  bathythermograph  data.    Magruder  [1970]  suggests  that 
Expendable  Bathythermograph  (XBT)  data  might  be  useful. 

The  purpose  of  this  thesis  is  to  see  if  information  about  the  distri- 
bution of  micros  true  ture  can  be  determined  from  XBT  traces.    If  the 
microstructure  distribution  can  be  determined  from  an  XBT  trace ,  then 
using  some  of  this  available  data,  can  the  distribution  be  related  to  the 
climatic  regions  of  the  ocean  for  use  as  a  rational  basis  for  assessing 
the  acoustic  variability  in  these  regions?    Finally,  what  operational 
techniques  can  be  developed  using  knowledge  of  the  microstructure  deter- 
mined from  the  XBT  trace  that  would  be  of  benefit  to  the  Navy  in  the 
operating  environment? 

Tens  of  thousands  of  XBT  traces  are  on  file  at  the  National  Oceano- 
graphic  Data  Center  (NODC) ,  with  the  number  increasing  each  day.     Under 
the  Fleet  Numerical  Weather  Central  (FNWC) ,  Cooperative  Oceanographic 
Observation  Program ,  XBT  traces  are  also  available  from  commercial 
carrier  trans-atlantic  and  trans-pacific  routes.     While  the  distribution 
is  by  no  means  uniform  in  space  or  time ,  it  is  the  only  source  of  micro- 
structure  data  covering  the  world  oceans. 


Osborn,  T.R.  ,  Report  of  the  Symposium  on  Turbulence  in  the  Ocean , 
Institute  of  Oceanography,  University  of  British  Columbia,  Vancouver, 
British  Columbia  ,  p.   20  ,  July  1970. 


II.    THEORY  OF  ACOUSTIC  PROPAGATION  IN  A 
RANDOMLY  INHOMOGENEOUS  MEDIUM 

Sheehy  [1950]  published  data  showing  that  the  magnitude  of  the 
acoustic  pressure  maximum  amplitude  fluctuations  (Coefficient  of 
Variation)  for  24  kHz  direct  path  transmissions  increased  proportional 
to  the  square  root  of  the  range.    However,  it  should  be  noted  from 
Figure  1  that  the  scatter  of  data  is  rather  large. 


100 


250    500    1000 
RANGE,  YARDS 


2  500  4000 


Figure  1.     Coefficients  of  Variation  of  the  amplitudes  of  direct 
signals  plotted  against  the  range  to  the  receiver.     Sound  frequency, 
24  kHz  [after  Sheehy,  1950]. 
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Bergmann  [1946]  using  the  application  of  ray  optics  to  investigate 
the  propagation  of  radiation  in  a  medium  with  random  inhomogeneities  , 
derived  an  approximate  formula  which  showed  that  the  signal  level  fluctua- 
tions increased  in  proportion  to  the  threerhalves  power  of  the  range.     He 
pointed  out  the  variance  with  Sheehy's  results ,  but  offered  no  explanation 
for  the  differences. 

Liebermann  [1951],  using  an  autocorrelation  function  calculated  from 
measured  horizontal  small-scale  temperature  inhomogeneities,  showed 
that  the  refractive  effects  of  the  horizontal  inhomogeneities  produce  a 
local  warping  (focusing  or  defccusing  of  sound)  that  causes  the  acoustic 
intensity  to  fluctuate.    He  found  that  refraction  should  only  occur  if  the 
frequency  is  high  enough  to  make  the  dimensions  of  the  inhomogeneities 
greatly  exceed  the  wave  length  of  the  radiation. 

Mintzer  [1953a]  derived  a  first  order  approximation  for  the  acoustic 
pressure  amplitude  variation  of  an  acoustic  pulse  traveling  in  a  medium 
with  small  random  variations  in  the  refractive  index  that  are  a  result  of 
random  temperature  variations.    Assuming  that  the  pressure  amplitude 
fluctuations  are  small  enough  to  allow  the  use  of  a  single  scattering 
approximation  and  that  the  refractive  index  at  a  point  only  changes  during 
the  time  between  pulses  and  not  while  the  pulse  is  passing  the  point,  he 
solved  the  Born  approximation  to  the  wave  equation  for  the  pressure 
amplitude  of  the  acoustic  pulse.    Then  using  the  Coefficient  of  Variation, 
V,  defined  as: 
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y2=  flpi2)  -  ICp>I  2 ,  where 

l<p>l  2 

^  /=  ensemble  average 
P    =  mean  presure  amplitude 

and  by  averaging  the  pressure  amplitude  over  the  pulse  length  instead  of 
using  the  maximum  value  of  the  pressure  amplitude  as  Sheehy  did  with  his 
experimental  data,    Mintzer's  theoretical  analysis  yields  the  following, 
for  ranges  large  compared  with  the  correlation  distance: 

V2=  2k0  a  r/N(p)dp     ;  where 
*o 

N(  p)=  autocorrelation  function  of  the  refraction 
index  variation  taken  along  the  line  between 
source  and  receiver 
kQ=         wave  number  of  the  sound 
r  =         range  from  source  to  receiver 
P=  rms  refractive  index  variation 

Using  a  Gaussian  correlation  function  the  coefficient  of  variation 
becomes: 

^       2     2      \     % 
V  =  (     r    k0  a  a  )     r      ,  where  a  ~  correlation 

length 
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This  theoretical  result  is  in  agreement  with  Sheehy's  experimental 
results  [Figure  1], 

Later  Mintzer  [1953b]  determined  over  what  range  the  first  assump- 
tion (use  of  a  single  scatterer  approximation)  was  valid.    He  found  that 
for  large  ranges  compared  with  the  acoustic  wave  length  and  the  inhomo- 
geneity  size ,  and  a  short  pulse  length  compared  to  the  time  required  for 
a  refractive  index  change: 

(1)    The  Coefficient  of  Variation  of  the  pressure 
amplitude  is: 

9  J-     2       2 

V*  =    7T  2  kQ     a     ar ,    ka  »  1 

2         h     3        2     2 

V    =   7t   k0      a     a"  r  ,  ka «  1 


(2)    The  single-scattering  assumption  used  is  valid 
for: 


k0r«(k0aa2)  2 


Mintzer  also  made  a  comparison  with  ray  theory  and  found  that  for  a 
frequency  high  enough  to  satisfy: 


rA0a2«  1 


the  intensity  coefficient  of  variation  computed  by  wave  theory  and  ray 
theory  are  the  same. 
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In  a  third  paper,  Mintzer  [1954]  considered  the  second  basic  assump- 
tion, that  the  refraction  index  varies  only  during  the  time  between  pulses. 
He  found  that  "if  the  pulse  length  is  much  less  than  the  time  constant  of 
the  refractive  index  correlation  function,  the  coefficient  of  variation  is 
independent  of  the  pulse  length,  and  the  pulse-correlation  function  is 
found  to  be  identical  with  the  time- dependent  part  of  the  refractive  index 
correlation  function.  " 

Skudrzyk  [  1957]  discussed  the  pressure  amplitude  fluctuations  caused 
by  water  inhomogeneities  considered  to  be  spherical  "patches"  that  are 
always  large  in  comparison  with  the  acoustic  wavelength.    He  showed  that 
in  the  vicinity  of  the  receiver  pressure  amplitude  fluctuations  are  caused 
by  focusing  or  dcf  ocusing  from  the  lens  action  of  the  patches ,  and  beyond 
a  certain  range  the  fluctuations  are  caused  only  by  an  interference  field 
which  results  from  deflected  rays  arriving  at  a  point  with  phase  shifts 
due  to  unequal  paths  of  travel.    He  also  pointed  out  that  the  interference 
component  of  scattering  is  independent  of  the  shape  of  the  patches ,  while 
the  focusing  component  can  be  significantly  affected  by  the  shape.     For 
example ,  a  spheroidal  patch  having  a  vertical  axis  less  than  the  horizontal 
axis  will  have  a  larger  focusing  component  than  a  spherical  patch.    At  low 
frequencies  (below  20  kHz)  interference  scattering  dominates ,  while  at 


Mintzer,  D.  ,  "Wave  Propagation  in  a  Randomly  Inhomogeneous 
Medium.      Ill,"    T.  Acous.  Soc.  Am.  ,  v.   26,  p.  186,  March  1954. 
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higher  frequencies  the  interference  component  may  be  so  small  in 
comparison  with  the  focusing  component  that  at  ranges  greater  than  the 
focusing  range  the  Coefficient  of  Variation  is  almost  constant.    The 
focusing  range  (  rQ  )  is  the  range  beyond  which  the  focusing  does  not 
change  and  additional  pressure  amplitude  fluctuations  are  caused  by 
interference  only.    Skudrzyk  found  the  focusing  range  to  be  a  function 
of  frequency  and  path  size. 


rQ  =  k  a     ,  wh 


ere 


k  =  acoustic  wave  number 


a  =  effective  radius  of  patch 


This  result  is  pictured  graphically  in  Figure  2.    If  the  sonar  frequency 

and  the  patch  size  are  known ,  then  the  focusing  range  can  be  determined. 

(For  example ,  a  24  kHz  sonar  has  a  focusing  range  of  200  meters ,  if  the 

patch  size  is  assumed  to  be  140  centimeters). 

a *• 

100  kc/tr 


50  kc/t 


20  kc/t 


10  kc/s 


5kcA 


2  kc/s 


I  kc/t 


10  20 


Figure  2.     Focusing  Range  (meters)  as  a  Functi    i  of  Frequency 
(kHz)  and  Patch  Size  (meters)  [Skudrzyk,  1957]. 
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Skudrzyk  shows  that  the  Coefficient  of  Variation  of  the  pressure 

amplitude  due  to  the  focusing  effect,  Vc,  is: 

2  o 

V  £  =    a      (z\      where  a  is  the  relative  root-mean- 
2       \ay       square  deviation  of  sound  velocity 


while  the  Coefficient  of  Variation  due  to  the  interference  effect  is: 

% 

Beyond  the  focusing  range ,  scattering  is  due  to  both  focusing  and  inter- 
ference so  the  Coefficient  of  Variation  becomes: 

V2     =vl+  v\ 

p  £  i 

and  for  spherical  patches : 


V2=    a2^    k2a  (r-0.44ro) 


which  is  Mintzer's  [1953b]  result  except  for  a  factor  of  1/2.    If  the 
patches  are  spheroidal  rather  than  spherical ,  then  the  focusing  may  be 
increased  significantly.    Whitmarsh,  Skudrzyk,  and  Urick  [1957]  found 
the  focusing  to  increase  by  a  factor 

4  6 

Vp  (Spheroidal)  x2      +        X4      +       6 

V     (Spherical)  16 

X  is  the  ratio  of  vertical  to  horizontal  axis  of  the 
patch 
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Figure  3  depicts  this  ratio  graphically.    The  interference  component  of 
scattering  is  virtually  independent  of  the  vertical  axis  dimension. 


K»r- 


Figure  3.    Increase  in  focusing  due  to  spheroidal  shape  of  a 
patch  as  a  function  of  the  ratio  of  the  lengths  of  the  vertical 
(b)  to  horizontal  (a)  axis  [Whitmarsh ,  Skudrzyk,  and  Urick, 
1957]. 


Since  published  experimental  results  pertained  only  to  one  frequency 
and  did  not  include  temperature  and  acoustical  fluctuation  measurements 
made  as  one  observation,  Whitmarsh,  Skudrzyk,  and  Urick  [1957]  conducted 
a  series  of  experiments  using  25kHz  and  60  kHz  transmissions  to  provide 
additional  data  to  test  the  existing  theories  for  scattering    of  acoustical 
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energy  in  the  sea.    They  used  a  gradiometer  bar  suspended  from  a  buoy 
(a  12  foot  long  bar  with  9  thermistor  beads)  to  record  temperature 
gradients  as  small  as  0.  001°C  and  to  record  the  temperature  fluctuations 
as  a  function  of  time.    They  point  out  that  all  the  existing  theories  for 
scattering  of  sound  in  an  inhomogeneous  medium  are  based  on  the  assump- 
tion of  a  spherical  patch  thermal  inhomogeneity  (one  such  patch  is  illustra- 
ted in  Figure  4).    The  theories  further  assume  that  the  patch  distribution 
in  the  medium  can  be  represented  by  an  exponential  or  Gaussian  correlation 
function. 


VViTT  a(r/o)3/2 
FOCUSSING    RANGE 


INTERFERENCE    RANGE 


Figure  4.  Illustration  of  the  various  types  of  acoustic  scattering 
due  to  a  thermal  patch  in  the  sea  [Whitmarsh ,  Skudrzyk ,  and  Urick , 
1957]. 


The  temperature  microstructure  measurements  revealed  that  the  assumed 
correlation  function  is  invalid  and  in  fact  the  results  she      the  temperature 
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fluctuations  increase  with  the  2/3  power  as  the  thermistor  spacing 
increases.    This  is  the  prediction  of  the  Kolmogorov  theory  for  homo- 
geneous turbulence ,  see  Figure  5.    However ,  even  replacing  the  correlation 
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Figure  5.     Horizontal  temperature  distributions  as  measured  in  the 
sea  [after  Whitmarsh ,  Skudrzyk,  and  Urick,  1957], 


function  by  a  new  function  representative  of  the  Kolmogorov  distribution, 
the  theories  for  sound  transmission  can  still  be  considered  vaHd. 
Whitmarsh ,  Skudrzyk ,  and  Urick  concluded  from  their  experiments  that 
the  present  theories  are  capable  of  predicting  the  approximate  variation 
in  the  pressure  ampHtude  to  be  expected  for  acoustical  transmissions  in 
an  inhomogeneous  medium. 
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Sagar  [1960]  analyzing  the  experimental  data  he  collected  off  New 
Zealand  suggests  that  sound  intensity  fluctuations  which  obey  the  Mintzer 

•y  2        2 

dependence  on  range  (V    cc    rf    a  a)  are  related  to  insolation  affects. 
His  data  is  interpreted  to  show  a    diurnal  cycle  of  growth  and  decay  of 
the  temperature  microstructure  in  the  surface  layers  of  the  sea  with 
the  structure  most  developed  by  late  afternoon  and  then  disappearing  over- 
night. He  also  found  and  described  the  water  to  have  a  "sandwich"  or 
layered  structure  rather  than  the  classical  "patch"  structure. 
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in.    RECENT  MEASUREMENTS  OF  MICROSTRUCTURE 

A.     REVIEW 

Stommel  and  Fedorov  [1967]  used  a  salinity,  temperature,  and  depth 
recorder  (STD)  to  take  measurements  of  the  fine  structure  (micro- 
structure)  in  the  main  thermocline  off  Timor  in  the  Arafura  Sea  and 
Mindanao  in  the  Philippine  Sea.    They  were  investigating  reports  of  fine 
structure  which  began  to  accumulate  with  the  advent  of  an  instrument 
like  the  STD,  which  continuously  records  temperature  and  salinity  as  a 
function  of  depth.    They  found  the  vertical    microstructure  in  the  main 
thermocline  to  consist  of  hundreds  of  layers  (laminae)  of  rather  homo- 
geneous temperature  from  2  to  40  meters  thick  separated  by  a  thin  layer 
(2  to  5  meters  thick)  of  high  temperature  gradient.    As  an  example ,  at 
one  station  a  mean  thermal  gradient  of  0.  040°C/meter  was  found,  while 
in  a  homogeneous  layer  the  gradient  was  greater  than  0.  3 5° C /meter.    By 
tracing  some  of  the  well  defined  laminae  in  a  horizontal  direction,  they 
found  the  laminae  extended  from  2  to  20  kilometers.    Though  the  STD 
resolution  was  only  2  meters ,  Stommel  and  Fedorov  felt  the  micro- 
structure  existed  on  a  smaller  scale  and  that  further  investigation  would 
be  of  interest. 

Cooper  and  Stommel  [1968]  made  microstructure  measurements  in  the 
main  thermocline  off  Bermuda  using  an  STD.    They  too  found  a  layered 
structure  with  the  homogeneous  layers  from  3  to  5  metei     thick  and  the 
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high  gradient  layers  between  10  and  15  meters  thick.    In  the  layers  of 
high  gradient,  temperature  changes  as  high  as  0.  50°C  were  found.    The 
horizontal  extent  of  the  layers  was  between  400  and  1000  meters.    This 
type  structure  was  present  out  to  a  distance  of  120  kilometers  from 
Bermuda.    Based  on  their  measurements ,  Cooper  and  Stommel  suggested 
that  most  of  the  Sargasso  Sea  exhibited  the  same  characteristics  in  the 
main  thermocline. 

Woods  [1968a]  conducted  an  extensive  study  of  the  microstructure 
in  the  summer  thermocline  of  the  Mediterranean  Sea.    A  temperature 
gradient  meter  was  used  to  obtain  measurements  of  temperature  gradient 
while  being  lowered  and  temperature  while  being  recovered.    The  instru- 
ment consisted  of  two  thermistors  mounted  50  centimeters  apart  on  a 
vertical  fiber  glass  rod.     Only  one  of  the  thermistors  was  used  for  the 
temperature  measurements.    The  instrument  had  an  0.01  C  accuracy  and 
a  depth  limitation  of  50  meters.    Woods  found  layers  with  a  moderate 
temperature  gradient  that  were  thick  in  comparison  with  alternate  layers 
(sheets)  of  high  temperature  gradient.     Figure  6  is  a  plot  of  his  results 
for  one  sounding.    High  gradient  sheets  can  be  seen  at  depths  of  10 ,  19 , 
24,  32,  41,  and  50  meters.    By  injecting  dye  from  suspended  packets  into 
the  different  layers ,  Woods  was  able  to  associate  laminar  flow  to  thin 
high  gradient  layers  and  a  weak  turbulent  flow  to  the  thick  homogeneous 
layers.     The  results  of  these  two  different  type  flows  is  a  strong  sheer 
across  the  high  gradient  sheets. 
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Figure  6.     Temperature  and  temperature- gradient  sounding 
for  0630  24  August  1966  in  the  Mediterranean  Sea  [Woods,  1968a]. 


The  next  year  Woods  [1968b]  conducted  another  investigation  of 
vertical  heat  flow  through  the  upper  ocean  in  the  Mediterranean  Sea  off 
Malta.    Again  the  sheets  in  the  thermocline  were  found  as  can  be  seen 
in  Figure  7. 
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Figure  7.    Temperature  and  temperature- gradient  soundings 
for  Maltese  waters,  14  September  1967  [Woods,  1968b], 


Lovett  [1968]  took  temperature  profiles  in  the  California  current 
off  Baja  California  using  a  sound,  velocity,  temperature,  and  pressure 
recorder  (SVTP).    He  too  found  a  definite  layered  structure  in  all  the 
temperature  profiles.    He  suggests  that  the  sharp  gradients  found  in 
the  step-like  structure  may  have  an  effect  on  acoustic  scattering, 
especially  so  at  high  frequencies. 
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Osbom  [1969]  using  a  free  fall  instrument  called  a  thermal 
gradiometer  made  microstructure  measurements  in  the  San  Diego  Trough , 
Florida  current  and  in  the  central  Pacific  between  San  Diego  and  Hawaii. 
He  found  regions  of  high  gradient  to  be  about  1  meter  thick  and  have  a 
gradient  around  10  times  the  mean  gradient.  Figure  8  is  a  record  of  the 
results  he  obtained  with  the  gradiometer.    By  tracing  a  region  of  high 
gradient  horizontally  he  found  they  extended  up  to  750  meters  if  the 
region  is  in  the  seasonal  thermocline ,  but  only  a  few  hundred  meters  if 
the  region  is  at  a  depth  greater  than  400  meters.    Osbom  suggests  that 
the  microstructure  is  on  a  scale  of  centimeter  order. 

Neshyba,  Neal,  and  Denner  [1969]  made  measurements  of  the  temper- 
ature mir^ostruciure  in  the  Arctic  Ocean  using  an  XBT  lowered  mechan- 
ically from  Fletcher's  Ice  Island  (T-3).    The  microstructure  found  had 
a  step-like  structure  with  layers  that  were  extremely  homogeneous  as 
can  be  seen  in  Figure  9.    The  transistion  layer  was  found  to  be  less  than 
20  centimeters  thick.   Figure  9b  shows  a  section  of  the  trace  recorded 
at  high  gain  (the  lip  at  the  end  of  a  transition  layer  is  due  to  recorder 
overshoot).    The  steps  from  220  to  300  meters  were  found  to  be  very 
stable.     Figure  10  is  a  time  series  taken  with  the  XBT  modified  to  fall 
at  7  ft/sec  instead  of  the  normal  20  ft/sec ,  and  demonstrates  the  stabil- 
ity.   Below  350  meters,  the  homogeneous  layers  are  larger,  but  still  are 
separated  by  steps  of  temperature  change.      This  is  also  apparent  from 
Figure  10.     It  should  be  noted  that  these  XBT  recordings  are  independent 
of  ship  motion  since  the  ice  island  is  a  very  stable  platf  <  rm, 
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Figure  8.     Corrected  data  from  the  upper  portion  of  the 
water  column  in  the  San  Diego  Trough  on  18  September  1968. 
The  curve  on  the  right  is  the  continuation  of  the  curve  on 
the  left      [Osborn,  1969]. 
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Figure  9.     Vertical  profile  of  temperature  under  T-3 
84  38'N,  128  21.6'W,  19  March  1969).     Shaded  section  of 
profile  (a)  shown  as  observed  in  profile  (b)  to  left  [Neshyba, 
Neal,  and  Denner,  1969]. 


30 


\  2145  (3/17/69)       0300  (3/18) 


2033(3/18) 


Figure  10.    Time  series  taken  from  T-3  between  2145Z 
17  March  and  2033Z  18  March  1969.    Temperature  increases 
from  right  to  left  [Neshyba,  Neal,  and  Denner,  1969]. 
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Simpson  [1971] ,  using  a  Bathysonde  CTD  system ,  examined  the  density 
and  thermal  microstructure  in  the  western  Irish  Sea.  He  found  vertical 
layers  up  to  10  meters  thick  and  in  cases ,  less  than  1  meter  thick.  The 
average  gradient  was  as  low  as  0.  002°C/meter  in  layers  up  to  10  meters 
thick  and  as  high  as  2°C/meter  in  the  layers  of  sharp  gradient.  Figure  11 
gives  three  temperature  profiles  of  the  region.  Horizontally  the  layers 
extended  up  to  a  distance  of  3  kilometers. 
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Figure  11.     Unretouched  Bathsonde  profiles  from  the  western 
Irish  Sea  of  temperature  and  conductivity  showing  pronounced 
microstructure  [Simpson,  1971]. 


B.        MICROSTRUCTURE  MECHANISMS 

The  exact  nature  of  the  processes  producing,  maintaining,  or  destroy- 
ing the  thermal  microstructure  still  are  unclear  because  of  the  complexity 
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of  the  problem.    Some  mechanisms  have  been  proposed  to  explain  the  re- 
sults of  observations  and  experiments. 

1.  Insolation 

Insolation  (radiation  from  the  sun  and  atmosphere)  results  in  the 
heating  of  the  water's  surface.    The  heat  is  transported  downward  mainly 
by  mixing  caused  by  wave  action.    The  wind-wave  action  tends  to  eliminate 
the  surface  layer  micros tructure.    However,  Sagar  [1960]  suggests  that 
after  sunrise  the  insolation  may  be  large  enough  to  dominate  over  wave 
action,  and  so  a  layered  structure  could  be  transported  down  the  thermal 
gradient.    The  depth  that  the  layers  reach  is  dependent  on  wave  action. 
Waves  are  required  to  produce  the  mixing,  but  if  they  are  too  strong, 
then  vigorous  mixing  will  eliminate  the  layered  structure.    Sagar  [I960] 
found  that  if  cloud  cover  existed  or  if  it  was  during  a  period  of  darkness , 
the  micros  tructure  was  absent. 

2.  Double  Diffusioi). 

This  mechanism  results  from  the  differences  in  the  molecular 
dif  fusivity  of  temperature  and  salinity  which  are  related  by  an  approxi- 
mate 100:1  ratio.     Due  to  the  differing  dif  fusivity  in  regions  where  a 
warm  saline  mass  of  water  overlies  a  cooler  mass ,  the  masses  become 
unstable,  forming  a  layered  structure.    Much  evidence  from  ocean  obser- 
vations [Cooper  and  Stommel,  1968,  and  Tait  and  Howe,  1968],  a  theor- 
etical study  [Turner ,  1968] ,  and  laboratory  experiments  [Turner  and 
Stommel,  1964,  Turner,  1967,  and  Stern,  1969]  support  this  mechanism. 
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Stern  [1969]  feels  that  a  turbulent  mixing  process  is  taking  place  in  thick 
layers  that  are  separated  by  thin  discontinuity  layers. 

3.  Internal  Waves 

Woods  [1968a]  observed  internal  waves  along  an  individual  micro- 
structure  sheet  to  be  in  long  coherent  trains.    If  these  waves  become 
unstable ,  they  would  break,  giving  up  their  energy  to  turbulent  processes 
causing  the  sheets  to  be  eliminated  by  mixing.    Orlanski  and  Bryan  [1969] 
have  proposed  that  internal  waves  breaking  could  also  establish  a  step  struc- 
ture by  leaving  a  "blob"  of  water  behind  that  would  spread  out  into  a  thin 
layer. 

4.  Current  Shear 

Current  shear  can  cause  vertical  turbulent  mixing  when  the 
Richardson  Number  is  less  than  0.  25  [Miles  and  Howard,  1963].    This 
turbulence  occurs  as  a  patch  which  is  smeared  out ,  leaving  a  scar  that 
can  lead  to  micros tructure. 

Lovett  [1969]  and  Pingre  [1969]  have  proposed  that  these  isolated 
turbulent  patches  can  also  be  spread  horizontally  by  turbulence  and  thus 
form  a  sheet.     Dye  studies  illustrate  this  layering  in  the  deep  ocean 

[Schuert,  1969]. 

The  complexity  of  the  problem  is  apparent  when  it  is  realized 
that  not  one ,  but  probably  a  combination  of  these  mechanisms  is  occurring 
simultaneously . 
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IV.    INFORMATION  CONCERNING  THE  THERMAL  MICROSTRUCTURE 
THAT  CAN  BE  GATHERED  FROM  EXISTING  EXPENDABLE 
BATHYTHERMOGRAPH  (XBT)  RECORDS 

A.  AVAILABLE  DATA 

Introduction  of  the  XBT  system  provided  a  method  relatively 
inexpensive ,  yet  at  the  same  time  accurate ,  of  obtaining  temperature- 
depth  profiles  originally  developed  for  ASW  operations,  the  system's 
application  has  been  expanded  to  include  Navy  research  ships  and  merchant 
ships  equipped  under  the  Fleet  Numerical  Weather  Central  (FNWC) , 
Cooperative  Oceanographic  Observation  Program.    The  commercial  ship 
routes  of  the  Pacific  and  Atlantic  Oceans  are  a  principal  source  of  ocean 
temperature  structure  data.    XBT's  taken  during  a  cruise  are  forwarded 
at  the  end  of  the  cruise  to  FNWC ,  Monterey  for  digitizing  and  then 
archived  for  further  use.     One  copy  of  the  digitized  record  along  with  the 
original  XBT  traces  is  sent  to  the  National  Oceanographic  Data  Center 
for  permanent  storage  in  their  archives.    As  of  May  1971,  approximately 
60  ,000  XBT's  had  been  digitized  by  FNWC  ,  Monterey  and  some  28  ,000  of 
these  were  from  the  northern  Pacific  area  [Dale,  1971]. 

B.  EXPENDABLE  BATHYTHERMOGRAPH  SYSTEM 

The  expendable  bathythermograph  system  currently  in  use  is  manu- 
factured by  the  Sippican  Corporation  of  Marion,  Massachusetts.    The 
system  is  composed  of  the  probe  (the  expendable  part  of  the  system) , 
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a  launcher,  and  a  recorder.    It  is  an  automatic  system  for  measuring 
and  recording  sea  temperature  as  a  function  of  depth.    The  specifications 
for  the  system  are  given  in  Table  I  [Sippican  Corporation,  1968]: 


TABLE    I 


Temperature  range 

Temperature  accuracy 

Depth  range 

Depth  accuracy 

Deployment  cycle  time 
Gradient  resolution 


Full  scale  stylus  travel 


28°F  to  96°F 
-2°C  to  35°C 

+  0.  4°F 
+  0.  2°C 

0  to  1500  feet 

0  to  2500  feet 

+  2%  or  15  feet,  whichever 
is  greater 

90  seconds  for  1500  ft  drop 

63%  of  a  step  change  in 
temperature  in  3  feet; 
95%  of  a  step  change  in 
9  feet 

1  second 


The  probe  is  released  from  a  cannister  which  pays  out  wire  to 
compensate  for  the  movement  of  the  launcher  platform  while  the  probe 
itself  pays  out  wire  to  compensate  for  its  fall.    The  probe  falls  with  a 
velocity  of  20  ft  /sec ,  and  has  been  found  to  have  an  average  terminal 
velocity  of  19.9  +0.7  ft/sec  [Gouzie,  Sanders,  and  Littlehale ,  1966]. 
A  thermistor  is  mounted  in  the  nose  of  the  probe ,  and  as  the  probe 
descends  through  the  water  column ,  the  water  temperature  is  sensed 
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as  it  flows  past  the  thermistor  and  out  the  tail  of  the  probe.    The 
response  time  for  the  thermistor  is  110  milliseconds  nominal  and  130 
milliseconds  maximum  [Choate ,  1970].      The  temperature  resolution 
of  the  probe  is  dependent  on  the  magnitude  of  the  temperature  change 
over  a  certain  depth  interval. 

The  recorder  uses  a  pressure  sensitive  paper  for  recording  the 
temperature  profile.    The  paper  moves  vertically  at  a  constant  rate 
proportional  to  the  fall  rate  of  the  probe.    The  chart  paper  displays 
311.656  ft/in  or  92.  307m/in  and  9.714°F/in  or  5.  286°C/in,  giving  a 
resolution  of  2  ft  and  0.1   F  for  visual  redding. 

C.    DIGITIZING 

The  digitizing  process  used  by  FNWC ,  Monterey  is  completely 
described  by  Bauer  [1969].     The  CALMA  Model  480  digitizer  is  a  X  -  Y 
recorder  that  records  at  least  every  0.15  inch  of  movement  in  either 
the  X  or  Y  direction.    The  speed  the  stylus  is  moved  down  the  trace 
determines  the  number  of  points  recorded.    The  maximum  resolution 
is  0.  01  inches  of  movement  in  the  X  or  Y  direction  which  is  equivalent 
to  .  053°C  or  .  097°F  in  temperature  and  0.  92  meters  or  3. 12  feet  in 
depth.    Each  trace  is  digitized  twice. 

The  digitized  data  is  further  processed  using  a  DIGITBT  computer 
program.     The  program  first  puts  the  data  in  a  common  format.     Then 
the  two  digitized  records  are  compared  point  by  point  in  the  X  and  Y 
dimensions.    If  the  maximum  differences  in  both  directions  are  less 
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than  0.035  inches  (0.186°C  or  0.  340°F  in  temperature  and  3.  22  meters 
or  10.  92  feet  in  depth) ,  the  two  records  are  averaged  on  a  point  by  point 
basis  to  produce  an  average  trace  which  is  used  for  the  remainder  of  the 
DIGITBT  program.    If  the  deviation  in  either  direction  is  greater  than 
0.  035  inches ,  then  the  records  are  rejected  and  must  be  redigitized. 

The  program  next  selects  points  that  will  be  printed  out  and  stored 
on  the  XBT  record  tapes.    These  points  are  called  INFLECTION  POINTS. 
Both  the  first  and  last  points  of  the  record  are  automatically  selected 
as  inflection  points.    The  other  points  of  the  record  are  screened  in 
"batches".     For  each  batch  the  observation  is  found  that  is  most  differ- 
ent from  that  expected  by  linear  interpolation  between  the  end  points 
of  the  batch.    If  that  difference  exceeds  0.015  inches  (0.079°C  or 
0. 146°F  and  1.  38  meters  or  4.  68  feet)  then  the  observation  is  selected 
for  output.    Otherwise  the  batch  expands  to  include  the  next  value  and 
the  inflection  test  is  repeated  until  an  inflection  is  found  and  until  the 
entire  record  has  been  screened.     Any  point  that  can  be  reconstructed 
by  linear  interpolation  is  not  required  and  therefore  is  removed  from 
the  output.     The  DIGITBT  program  converts  the  inflection  points  from 
inches  to  meters  and  Centigrade  degrees  for  output.    The  depth  is 
rounded  to  the    nearest  whole  unit  and  temperature  to  the  nearest  tenth 
of  a  degree.     Figure  12  is  an  example  output  from  the  DIGITBT  program 
and  is  the  information  stored  in  FNWC,  Monterey's  records.    The  title 
line  contains  a  ship  identification  code  ,  day ,  month  ,  year ,  time , 
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Figure  12.     Digitized  record  [Dale  and  Stevens,  1970  ]. 

latitude ,  longitude  ,  base  line  temperature  ,  grid  type ,  and  maximum 
depth  of  probe  (15  for  1500  ft  XBT  or  25  for  2500  ft  XBT).     Figure  13 
is  a  comparison  of  the  original  XBT  trace  and  the  reconstructed  trace 
from,  a  DIGITBT  record  similar  to  the  example  in  Figure  12.    As  could 
be  expected,  the  digitizing  process  has  eliminated  some  of  the  small 
scale  micros  true  ture. 

For  this  thesis  the  data  used  is  from  the  FNWC,  Monterey's  data 
bank  that  has  been  processed  using  the  DIGITBT  program. 
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Figure  13.    Comparison  of  an  original  XBT  trace  with  the  trace 
reconstructed  from  the  digitized  data.     The  digitized  trace  is  plotted 
to  a  scale  of  90  meters/in  due  to  computer  requirements  while  the 
original  trace  is  scaled  to  92.  3  meters/in. 
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D.      THERMAL  MICROSTRUCTURE 

Thermal  microstructure  for  the  purpose  of  this  thesis  is  defined 
as  the  temperature  fluctuations  as  a  function  of  depth  that  exist  on  a 
scale  that  is  smaller  than  the  scale  of  the  gross  temperature  structure 
■which  is  derived  from  the  large  scale  processes  such  as  weather,  latitude, 
season,  and  water  mass  transport. 

To  analyze  the  XBT  trace  for  microstructure ,  a  method  was  sought 
to  separate  the  temperature  perturbations  (microstructure)  from  the 
mean  of  the  XBT  trace.    A  method  to  filter  the  trace  and  subtract  the 
mean  portion  of  the  curve  was  desired.    Several  type  filters  were  consid- 
ered.   A  binomial  filter  was  selected  primarily  because  Roden  [1971]  had 
applied  a  modified  binomial  filter  to  similar  data  with  favorable  results  = 

Roden  utilized  a  low  pass  filter  to  develop  a  mean  series.    This  mean 
series  was  then  subtracted  from  the  original  series  to  give  a  temperature 
perturbation  series.    The  infrmation  was  therefore  effectively  passed 
through  a  high  pass  filter.    Using  this  scheme  the  temperature  pertur- 
bations are  considered  to  be  associated  with  the  thermal  microstructure. 
If  the  original  series  is  denoted  by  s(z),  then  the  mean  series,  s*(z),  is 

found  by: 

m 

2 


s*  (z)  =  2_        w(i)s(z+1) 


m 
i=-T 
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where  w(i)  is  a  smoothing  array  that  has  the  following  properties: 

w(i)  =  w(-i) 

£>(£)  =  1 

i 

The  smoothing  array  is  found  by  evaluating  the  below  modified  binomial 
sequence. 


w(i)  = 


m! 


(m-i)!i! 


The  parameter,  m,  is  even  and  controls  the  width  (number  of  weighted 
points)  of  the  array.    The  parameter,  p,  is  an  arbitrary  exponent  which 
controls  the  peakedness  of  the  array.    Roden  found  that  a  statistically 
stable  mean  series  could  be  obtained  if  m  -  40 ,  with  p  =  0.5.    Increasing 
m  produces  more  smoothing,  but  the  series  is  physically  shortened  as  will 
be  seen  later.     The  response  of  the  modified  filter,  r(k) ,  is  given  by 

r(k)  =  l-/Jw(i)cos(ik  A  ft  ) 
i 

where  k  =  wave  number  (    Vwave  length)  and     A#  =  filtering  interval 

(distance  between  the  41  weighted  points).     Figure  14  is  the  solution  for 

the  response  function  for  filtering  intervals  (  A  fi  )  of  1 ,  2  ,  and  3  meters. 

The  smaller   A  J£  the  sharper  is  the  response  curve  at  a  cut-off  wave 

length,  however,  the  cut-off  wave  length  is  lowered.    It  should  be 

remembered  that  with  a  binomial  filter  the  amplitude  response  is  not  very 

predictable  in  the  region  of  cut-off. 
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Figure  14.     Computer  solution  for  modified  binomial  filter 
with  filtering  intervals  of  1 ,  2 ,  and  3  meters. 
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A  theoretical  XBT  trace  was  constructed  for  test  purposes.    An 
isothermal  layer  of  50  meters  was  followed  by  a  thermocline  simulated 
with  an  exponential  decay  function  (e        )  that  decayed  in  650  meters. 
The  microstructure  was  constructed  by  superimposing  a  sine  wave 
function  (sin2  w  kz)  on  the  exponential  function.    The  amplitude  was 
varied  to  represent  the  temperature  fluctuations  desired,  and  the  wave 
length  varied  to  represent  the  microstructure  layer  thickness  scale. 
Figure  15  is  a  composite  picture  of  a  theoretical  trace,  the  mean 
filtered  trace ,  and  the  resulting  temperature  perturbation  curve  for  a 
7  meter  wave  length,  0.  5°C  temperature  variation,  and  a  filtering 
interval  of  1  meter.    As  expected,  the  perturbation  curve  is  the  input 
sine  wave  curve  with  both  the  wave  length  and  amplitude  of  the  temper- 
ature variation  being  preserved.     Figure  16  is  another  set  of  theoretical 
curves  showing  that  for  the  one  meter  filtering  interval  if  the  wave 
length  is  34  meters  and  the  temperature  amplitude  variation  is  again 
0.  5°C  the  temperature  perturbation  amplitude  has  been  amplified,  while 
the  wave  length  is  preserved.    It  should  be  noted  that  the  34  meter  wave 
length  is  on  the  slope  part  of  the  response  curve  in  Figure  14 ,  and  this 
shows  the  instability  of  the  filter  in  this  region  as  predicted.     If  the 
filtering  interval  is  increased  to  3  meters  then  the  original  sine  wave 
amplitude  will  again  be  passed. 


44 


TEMPERATURE 


CO 


100 


200 


XBT 
10  15 


20 


PERT 
10  0 


10 


__^- 


«-^ 

300 

U) 

« 

LLI 

H- 

UJ 

2 

^^ 

400 

X 

►- 

o. 

Ul 

Q 

500 

-=5£ 


600 


700 


800 


:XBT 


MEAN 


;-     vr  -r^  peri 


. 


10  15 

MEAN 


20 


Figure  15.     Result  of  modified  binomial  filter  applied  to  Theorectical 
XBT  trace.     From  left  to  right,  Theoretical  XBT  trace,  Mean  curve, 
and  Temperature  Perturbation  curve  for  a  7  meter  wave  length,  0.  5°C 
Temperature  variation,  and  a  1  meter  filtering  interval. 
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Figure  16.     Result  of  modified  binomial  filter  applied  to  theoretical 
XBT  trace  showing  the  amplification  in  the  temperature  perturbation 
amplitude  when  the  wave  length  of  the  perturbation  is  on  the  slope 
part  of  the  response  curve  of  Figure  14.     From  left  to  right,  Theo- 
retical XBT  trace,  Mean  curve,  and  Temperature  Perturbation  curve 
for  a  34  meter  wave  length,  0.  5°C  Temperature  variation,  and  a 
1  meter  filtering  interval. 
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E.      DATA  SELECTED  FOR  ANALYSIS 

A  selection  of  129  XBT  traces  has  been  analyzed  for  this  thesis  using 
Roden's  modified  binomial  filter  with  1,2,  and  3  meter  digitizing  inter- 
vals.    All  of  the  traces  were  from  the  northern  Pacific  Ocean.     Digitized 
information  was  furnished  by  the  Fleet  Numerical  Weather  Central  (FNWC) , 
Monterey ,  and  was  processed  using  the  IBM  360/67  at  the  Naval  Post- 
graduate School's  W.   R.   Church  Computer  Center.     The  traces  were 
extracted  from  the  FNWC  data  bank  by  selecting  geographic  grid  squares. 
The  selection  of  these  grid  squares  was  based  on  three  criteria.    The  first 
and  most  important  was  to  get  a  sampling  from  the  Subtropic ,  Subarctic , 
Kuroshio  Transport  and  Region  of  Confluence  climatic  regions  proposed  by 
Tully  [1964].     Secondly  it  was  desired  to  secure  data  for  comparison  with 
copies  of  original  XBT  traces  obtained  from.  J.  F.  T.  Saur  of  NOAA's 
National  Marine  Fisheries  Service  at  the  Fishery- Oceanography  Center  in 
La  Jolla,  California.    The  third  criteria  was  to  request  a  reasonable  quan- 
tity of  data  from  FNWC.     Data  was  desired  for  all  months  of  the  year  in 
each  grid  square  ,  but  unfortunately  the  coverage  fell  short  of  expectations. 
Another  shortcoming  in  the  data  is  that  the  Subarctic  region  has  not  been 
sampled  as  often  with  XBT  as  the  Subtropic  region  and  therefore ,  data  is 
not  very  plentiful  from  this  area. 
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F.      DESCRIPTION  OF  MICROSTRUCTURE  ON  XBT  TRACE 

Figure  17  is  a  composite  diagram  of  an  actual  XBT  trace,  the 
digitized  trace,  a  mean  curve  found  with  Roden's  filter,  and  the  temper- 
ature perturbation  curve  using  a  1  meter  filtering  interval.    The  scale  of 
the  actual  XBT  trace  is  92.  3  meters/in  while  the  other  three  traces  are 
scaled  90  meters/in  for  computer  requirements.    By  comparing  the  actual 
XBT  trace  and  the  digitized  reproduction  it  can  be  seen  that  some  of  the 
very  fine  scale  temperature  micros  true  ture  has  been  eliminated  by  the 
digitizing  process  and  some  larger  scales  appear  to  have  been  amplified 
by  the  process.    The  exactness  of  the  digitized  trace  is  considered  to  be 
sufficiently  reliable  to  achieve  results  that  are  realistic  and  comparable. 

The  first  and  last  20  meters  of  the  mean  trace  are  lost  because  the 
filter  used  20  points  at  1  meter  intervals  on  both  sides  of  a  depth  to 
compute  the  mean  value  for  that  depth.    If  a  2  meter  filtering  interval 
is  used,  then  a  total  of  80  meters  is  lost  and  the  3  meter  interval  reduces 
the  length  by  120  meters.    Since  the  XBT  record  is  comparatively  short 
to  begin  with,  a  750  meter  trace  would  be  reduced  by  16  percent  when  using 
the  3  meter  interval,  and  if  the  trace  is  a  1500  ft  XBT  (460  meters)  26 
percent  of  the  record  will  be  lost. 

By  closely  comparing  the  digitized  XBT  trace  and  the  temperature 
perturbation  curve  in  Figure  17 ,  it  can  be  seen  that  any  time  an  inflection 
point  (as  defined  in  Section  IV  C)  is  encountered  on  the  XBT  trace ,  a  change 
in  the  slope  results  on  the  temperature  perturbation  curve.     The  number  of 
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Figure  17.     Composite  diagram  of  the  original  XBT  trace,  the 
reconstructed  Digitized  trace ,  the  Mean  curve  found  by  applying 
Roden's  filter,  and  the  Temperature  Perturbation  curve.    The 
original  trace  is  scaled  92.  3  meters/in  while  the  remaining  curves 
are  scaled  90.  0  meters/in. 
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slope  changes  can  be  counted  and  divided  into  the  length  of  the  trace 
producing  a  mean  temperature  perturbation  vertical  thickness  scale  which 
is  representative  of  the  microstructure  mean  vertical  scale.    The  same 
value  can  be  obtained  by  dividing  by  the  number  of  inflection  points  that 
was  selected  for  the  digitized  record.    To  be  sure  that  only  the  micro- 
structure  is  being  considered,  the  inflection  points  or  slope  changes 
above  the  seasonal  thermocline  must  not  be  used  nor  the  depth  interval 
from  the  surface  to  the  seasonal  thermocline.     This  must  be  determined 
for  each  trace  independently.     The  amplitude  of  the  temperature  pertur- 
bations is  defined  as  the  magnitude  of  the  temperature  difference  between 
the  mean  curve  and  the  XBT  trace. 

Fiaure  18  demonstrates  how  the  temperature  perturbations  are 
affected  by  changing  the  filtering  interval.     Examination  of  the  three 
traces  below  100  meters  shows  that  if  the  filtering  interval  is  increased, 
then  the  shape  of  the  perturbation  curve  is  changed.    In  general ,  it  was 
found  that  the  mean  temperature  perturbation  vertical  thickness  scale 
was  increased  when  the  filtering  interval  (  A$)  was  increased,  and  the 
temperature  amplitude  was  increased,  but  not  at  all  depths.    It  should 
also  be  noted  that  more  of  the  trace  is  influenced  by  the  mixed  layer  depth 
as  the  interval  is  increased,  so  that  for  the  3  meter  curve  150  meters 
are  not  useable  in  computing  the  mean  temperature  perturbation  vertical 
thickness  scale. 
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Figure  18.     Temperature  perturbation  curves  for  a  filtering 
interval  of  1 ,  2  ,  and  3  meters. 
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Previous  observations  of  the  microstructure  indicate     scales  below  a 
few  meters  are  important.    Since  increasing  the  filtering  interval  pro- 
duced an  increase  in  the  mean  temperature  perturbation  vertical  thickness 
scale ,  the  perturbation  amplitude  and  the  variance ,  it  was  decided  that 
the  1  meter  filtering  interval  was  more  effective  for  processing  the 
digitized  XBT  data. 

G.     OCEANOGRAPHIC  REGIONS  OF  THE  NORTH  PACIFIC  OCEAN 
1.      Description 

As  previously  mentioned,  Tully  [1964]  divided  the  north  Pacific 
Ocean  into  seven  regions  to  be  used  for  examining  and  forecasting  oceano- 
graphic  characteristics  and  behavior  on  a  similar  scale  as  the  Weather 
Services  have  been  doing  for  many  years.    In  defining  the  regions ,  he 
considered  volume  transport,  heating  and  cooling,  precipitation  and  evap- 
oration, mixing  (wind  and  convection),  internal  waves,  and  the  effects 
these  processes  have  on  changing  the  properties  of  sea  water.     Figure  19 
is  a  diagram  of  the  Pacific  Ocean  and  the  seven  regions.    The  Tropic, 
Subtropic,  Arctic,  and  Subarctic  regions  are  the  water  mass  formation 
areas  where  due  to  slow  and/or  zonal  transport ,  the  water  mass  has  time 
to  establish  the  regional  properties.     The  normal  time  required  to  achieve 
these  properties  is  one  heating  and  cooling  season. 

The  Tropic  region  is  bounded  by  Latitudes  10°S  and  10°N.    Wind  is 
the  dominant  mixing  factor.    The  temperature  is  warm  with  a  small  season- 
al variation.     The  flow  is  generally  zonal  with  a  net  westwr-d  transport. 
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Figure  19„     Regions  of  the  North  Pacific  Ocean  [Tully,  1964]. 
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Figure  20.     Drift  in  the  North  Pacific  Ocean  [Tully     1964]. 
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The  Subtropic  is  the  largest  of  the  regions.    The  flow  is  zonal 
and  relatively  slow,  as  is  shown  in  Figure  20.    The  evaporation  exceeds 
precipitation,  causing  an  evaporation- driven  convection  mixing.    The  density 
structure  is  primarily  a  function  of  the  temperature  structure.     From 
spring  to  autumn  more  heat  is  added  during  the  day  than  is  lost  during 
cooling  at  night,  so  that  in  September  the  warmest  surface  water  occurs. 
The  region  extends  north  to  between  Latitudes  40°N  and  45°N. 

The  Subarctic  region  is  characterized  by  excess  precipitation 
and  as  a  result  there  is  a  distinctive  salinity  structure.     Fresh  water 
accumulates  in  the  upper  100  meters  with  a  halocline  in  the  next  100 
meters  where  salinity  increases  about  1  percent.     The  density  in  this 
region  is  primarily  a  function  of  salinity.    Mixing  is  due  to  the  wind  alone. 

The  Pacific  Arctic  region  also  has  an  excess  of  precipitation 
and  a  salinity  halocline.    There  is  an  ice  cover  in  the  winter  and  until  mid- 
summer the  heat  added  is  used  to  melt  the  ice ,  followed  by  a  rapid  and 
shallow  warming  due  to  the  long  northern  days.     Convective  mixing  is 
dominant  in  this  region. 

The  Transport  regions  occur  at  the  western  boundaries  of  the 
ocean.    An  intensified  meridianal  flow  transports  water  between  the 
regions  of  zonal  flow  (Figure  20).    The  water  does  not  remain  in  these 
transport  regions  long  enough  to  fully  acquire  the  regional  characteristics. 

The  Kuroshio  Current  is  the  western  boundary  for  the  Sub- 
tropic  region.    It  transports  a  thick  layer  of  warm  saline  water  formed 
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in  the  Tropics  near  the  Philippine  Islands  past  the  coast  o£  Japan  losing 
some  of  its  heat  along  the  way ,  but  never  becoming  fully  adjusted  to  the 
climate  of  the  region.    Its  boundaries  are  easily  identified  due  to  large 
horizontal  temperature  gradients. 

The  Oyashio  Current  is  the  western  boundary  for  the  Sub- 
arctic region.    It  has  a  smaller  rate  and  volume  transport  than  the 
Kuroshio  Current.     The  flow  is  coast  hugging  and  becomes  strongest 
during  the  summer  when  it  is  swelled  because  of  the  ice  melting  in  the 
Arctic. 

The  Confluence  Region  is  the  area  where  the  Kuroshio  and 
Oyashio  Currents  meet.    The  resultant  flow  is  eastward.     The  boundaries 
between  these  two  flows  can  have  temperature  gradients  as  high  as  5°C 
in  a  mile.     Eddies  from  one  stream  can  be  dissolved  in  the  other  stream. 
As  the  Confluence  moves  east  it  looses  its  identity  and  is  dissipated  in 
the  slow  eastward  trans-ocean  flow  (Figure  20).     More  than  two  years 
are  required  for  the  water  to  cross  from  Asia  to  America  [Uda,  1963] 
so  the  surface  processes  can  fully  re-establish  the  regional  characteristics. 
2.      Microstructure  Distribution 

The  data  available  to  this  author  necessitated  that  the  regions 
investigated  be  restricted  to  the  Subtropic  ,  Subarctic  ,  Kuroshio  Current 
and  Confluence  regions.     The  first  region  and  the  one  most  thoroughly 
examined  was  the  Subtropic  for  the  purpose  of  comparing  the  results  with 


55 


those  of  Roden  [1971]  who  had  conducted  a  detailed  study  of  the  vertical 
temperature  and  salinity  structure  in  the  central  North  Pacific  in  April 
1968  using  STD  data. 

Six  locations  were  selected  in  the  Subtropic  region  to  see  if 
the  mean  temperature  perturbation  vertical  thickness  scale  was  related 
within  the  region.     The  mean  temperature  perturbation  vertical  thickness 
scale  was  found  for  each  of  the  68  available  traces  using  a  1  meter  filter- 
ing interval.    The  average  was  then  found  for  the  traces  from  each  of  the 
six  locations.    This  average  mean  temperature  perturbation  vertical 
thickness  scale  will  be  referred  to  as  the  average  thickness  scale.    Also 
for  each  trace ,  the  maximum  amplitude  of  the  temperature  perturbation 
was  found  and  then  averaged  for  each  location.    Table  II  presents  the 
results  for  the  Subtropic  region.    The  minimum,  maximum,  and  average 
values  for  each  location  are  listed  along  with  the  number  of  traces  analyzed 
for  each  location.     Roden' s  values  are  given  for  similar  locations  for 
comparison.    The  average  thickness  scale  varies  between  8.7  and  15.0 
meters  which  compares  favorably  with  Roden' s  results  which  show  a  varia- 
tion from  13  to  16  meters  for  this  region.    One  reason  that  a  wide  range 
of  variations  was  found  is  that  in  analyzing  the  data ,  Roden  used  zero 
crossings  to  compute  the  thickness  scale  while  in  this  analysis ,  inflection 
points  from  the  digitized  record  were  used  because ,  as  was  pointed  out  in 
Section  IV  F,  they  represent  layers  resolved  by  the  digitizing  process.     In 
addition ,  the  68  traces  that  were  averaged  represent  all  of  the  months  of 
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the  year  and  not  just  April.    Another  factor  is  that  the  maximum  depth 
of  the  XBT  traces  was  750  meters  while  the  STD  record  extended  to 
1500  meters. 

The  average  maximum  amplitude  of  the  temperature  perturba- 
tions varies  between  0.  30  to  0.  49°C  which  is  from  the  same  minimum  to 
three  times  the  maximum  value  Roden  found.    This  difference  is  probably 
largely  due  to  the  STD  not  recording  the  total  magnitude  of  a  temperature 
change  since  its  response  time  is  about  one-third  as  fast  as  the  XBT. 
The  STD  is  also  temperature  compensated  in  the  salinity  sensor  and  large 
temperature  changes  in  the  water  column  cause  a  salinity  spike  that  will 
show  by  comparison  with  the  temperature  record  the  whole  change  was 
net  recorded.     Another  important  factor  a_f:f ecting  the  STD  is  that  it  is 
not  free  from  the  vertical  motion  of  the  platform.     This  vertical  motion 
as  well  as  the  physical  size  of  the  instrument  could  distort  the  temper- 
ature structure  before  it  is  sampled. 

Table  III  gives  the  results  for  the  Subarctic  region.     Only  12 
traces  were  analyzed.     The  presentation  is  the  same  as  for  the  Subtropic 
region  with  the  exception  of  the  omission  of  the  Roden  values  since  he  did 
not  have  any  results  from  this  region.     The  average  thickness  scale  varied 
between  17.  4  and  82  meters.    The  location  for  the  maximum  thickness 
scales  were  in  the  Gulf  of  Alaska  in  the  Alaskian  gyre.    The  three  traces 
analyzed  were  for  March,  April,  and  July  with  a  mean  thickness  scale  of 
84  m  ,    101m,  and  51m  respectively.    The  general  behavior  (more 
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micros  true  ture  activity  present  during  the  long  summer  days  and  less 
during  March ,  the  coldest  month)  of  the  conditions  is  as  would  be  expected 
since  the  area  is  unstable  and  with  the  cold  winter,  convection  mixing  would 
tend  to  homogenize  the  water  column  and  with  the  summer  heating  the 
micros  true  ture  would  reform.    The  average  maximum  amplitude  of  the 
temperature  perturbations  varied  between  0.16  and  0.  42°C. 

The  Kuroshio  Current  and  Region  of  Confluence  results  are 
shown  in  Table  IV.     The  only  comment  on  the  data  is  that  at  40 °N  149°F 
the  average  thickness  scale  is  smaller  than  for  the  rest  of  the  region. 
This  small  value  indicates  a  large  concentration  of  micros tructure  activity 
and  should  be  expected  because  this  is  the  area  where  the  Kuroshio  and 
Oyashio  currents  meet,  creating  a  region  of  intense  vertical  shear  that  is 
effective  in  producing  microstructure. 

Table •  V  gives  results  from  two  other  areas.    One  location  was 
in  the  Luzon  Strait  between  Luzon  and  Taiwan.     This  also  demonstrates 
that  there  is  a  large  amount  of  microstructure  present  and  again  is  in  a 
location  that  current  shear  could  be  expected  due  to  the  restriction  of 
flow. 

Also  in  Table  V  are  the  results  of  a  28  trace  time  series  taken 
on  two  days  in  the  California  Current  200  nautical  miles  west  of  Monterey , 
California.     Even  though  the  extremes  of  the  mean  thickness  scale  are 
widely  separated,  19  of  the  traces  or  68  percent  were  within  +  2  meters 
of  the  average  value  (13.0  meters).     Table  VI  gives  the  lo-al  time,  mean 
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thickness  scale ,  and  maximum  amplitude  of  the  temperature  perturbations 
for  each  of  the  28  traces.    The  hours  are  not  consecutive  because  some  of 
the  traces  were  rejected  by  FNWC,  Monterey  when  received  for  process- 
ing.   Without  weather  information,  it  is  not  possible  to  correlate  the 
microstructure  mean  thickness  scales  and  the  local  time. 

The  average  thickness  scale  data  from  Tables  II,  III,  IV,  and  V 
has  been  plotted  in  Figure  21  with  the  latitude  of  the  location ,  and  shows 
quite  clearly  that  in  the  north  Pacific  Ocean  for  an  increase  in  latitude 
the  average  mean  temperature  perturbation  vertical  thickness  scale  in- 
creases.   A  similar  plot,  Figure  22,  shows  that  the  average  maximum 
amplitude  of  the  temperature  perturbations  decreases  with  an  increase  in 
latitude.     Due  to  the  scatter  of  data  the  relationship  is  net  as  clear  as 
for  the  thickness  scale.     Figure  23  is  a  plot  of  the  maximum  amplitude 
of  the  temperature  perturbations  versus  the  average  mean  temperature 
perturbation  vertical  thickness  scale.    Although  the  scatter  of  data  is 
wide ,  the  relationship  of  the  plot  can  be  determined.     For  an  increase  in 
the  thickness  scale  of  the  microstructure  the  maximum  amplitude  of  the 
temperature  perturbation  will  decrease.     Sverdrup ,  Johnson,  and  Fleming 
[1942]  have  shown  an  increase  in  water  instability  when  preceding  north 
from  Hawaii  to  Alaska.    With  an  increase  in  instability  the  thick  layers 
of  turbulent  convective  mixing  and  of  small  temperature  gradient  would 
become  thicker,  eliminating  the  thin  layers  of  high  gradient.     This  is  a 
possible  explanation  for  the  results  shown  in  Figure  21. 
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Climatological  data  was  not  obtained  so  that  no  conclusions 
could  be  drawn  concerning  the  local  time  and  the  microstructure  activity. 
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TABLE  V 


Other  Areas 


California  Current 

36-37N, 

125-126W 

28  traces 

Mean  Temperature 

MIN 

7.9 

Perturbation  Vertical 

MAX 

19.3 

Thickness  Scale  (m) 

AVERAGE 

13.0 

Maximum  Amplitude 

MIN 

0.16 

Of  The  Temperature 

MAX 

0.33 

Perturbations  (°C) 

AVERAGE 

0.24 

Luzon  Strait 

20-24N, 

120-124E 

2  traces 

Mean  Temperature 

MIN 

6.5 

Perturbation  Vertical 

MAX 

7.8 

Thickness  Scale  (m) 

AVERAGE 

7.15 

Maximum  Amplitude 

MIN 

0.43 

Of  The  Temperature 

MAX 

1.00 

Perturbations  (°C) 

AVERAGE 

0.71 
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TABLE  VI 


California  Current  Time  Series 


Local 

Mean  Temperature 

Maximum  Amplitude 

Date 

Time 

Perturbation  Vertical 

Of  The  Temperature 

Thickness  Scale  (m) 

Perturbations  (°C) 

11  Nov  68 

0200 

15.0 

0.24 

0400 

15.4 

0.22 

0500 

12.8 

0.32 

0600 

13.3 

0.16 

0700 

16.6 

0.20 

0900 

13.1 

0.20 

1000 

15.9 

0.18 

1100 

15.0 

0.33 

1200 

10.7 

0.29 

1300 

11.0 

0.25 

13  Nov  68 

2300 

7.9 

0.20 

14  Nov  68 

0100 

12.5 

0.22 

0200 

8.7 

0.16 

0400 

19.3 

0.10 

0500 

13.5 

0.16 

0600 

15.7 

0.18 

0800 

10.7 

0.22 

0900 

13.1 

0.20 

1000 

9.9 

0.24 

1100 

13.5 

0.16 

1200 

12.9 

0.14 

1300 

11.0 

0.20 

1400 

10.5 

0.27 

1500 

13.3 

0.20 

1700 

13.1 

0.22 

1800 

14.3 

0.22 

1900 

12.8 

0.26 

2100 

12.4 

0.24 
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5  10  15  20  25 

AVERAGE  MEAN  TEMPERATURE  PERTURBATION  VERTICAL 
THICKNESS  SCALE  (meters) 


Figure  21.     Plot  of  the  Average  Mean  Temperature  Perturbation 
Vertical  Thickness  Scale  versus  Latitude  showing  an  increase  in 
the  vertical  thickness  scale  for  an  increase  in  latitude  between 
20°N  and  55°N.     (  A  ,  points  along  Longitude  156-15^W;  •  ,  all  other 
points). 
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Q2  04  06  0.8 

AVERAGE  MAXIMUM  AMPLITUDE  OF  THE  TEMPERATURE 
PERTURBATIONS  (°C) 


Figure  22.     Plot  of  the  Average  Maximum  Amplitude  of  the 
Temperature  Perturbations  versus  Latitude  showing  the  Amplitude 
decreases  with  an  increase  in  latitude  between  20°N  and  55  N. 
(  A  ,  points  along  Longitude  156-157°W;  •  ,  all  other  points). 
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AVERAGE  MAXIMUM  AMPLITUDE   OF   THE   TEMPERATURE 
PERTURBATIONS    (°C) 


Figure  23.     Plot  showing  the  Average  Maximum  Amplitude  of  the 
Temperature  Perturbations  decreases  as  the  Average  Mean  Temper- 
ature Perturbation  Vertical  Thickness  Scale  increases  between 
20°N  and  55°N.      (  A  ,  points  along  Longitude  156-157°W;  •  ,  all  other 
points). 
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V.    OPERATIONAL  ANALYSIS 

A.     OPERATIONAL  REQUIREMENT 

Underwater  acoustic  energy  transmissions  are  extremely  vital  to  the 
Navy,  especially  in  the  anti-submarine  warfare  (ASW)  field.    Some  of  the 
military  uses  and  some  of  the  non-military  uses  of  underwater  sound  as 
reported  by  the  National  Academy  of  Sciences  [1970]  are  listed  below  by 
frequency  bands. 

1.  Less  than  100  Hz 

a.  Passive  and  active  long  range  sonar 

b.  SOFAR  navigation  and  positioning 

c.  Deep  exploration  and  mapping 

2.  100  to  5000  Hz 

a.  Shallow  seismic  exploration  and  mapping 

b.  Warning  bells 

c.  Scuba  diver  communications 

d.  Passive  and  active  ship  sonars 

e.  Mines 

3.  5  to  15  kHz 

a.  Underwater  telephone  transmissions 

b.  Depth  sounding  equipment 

c.  Fish  finding 

d.  High  frequency  ship  sonar 
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4.  15  to  100  kHz 

a.  Short  range  underwater  telephone  transmissions 

b.  Depth  soundings 

c.  Fish  finding 

d.  Ship  positioning 

e.  Object  location 

f .  Low  resolution  scanning  sonars 

g.  High  resolution  sound  velocity  profiling 
h.      Homing  torpedoes 

5.  Greater  than  100  kHz 

a.  Doppler  sonar  for  navigation  or  docking 

b.  High  resolution  scanning  sonar 

For  low  frequency  acoustic  transmissions  thermal  microstructure 
should  not  have  a  significant  effect  on  the  propagation  because  the  acoustic 
wave  length  is  large  relative  to  the  effective  radius  of  a  patch  of  thermal 
inhomogeneities.     However,  losses  at  low  frequencies  have  been  related  to 
scattering  by  the  index  of  refraction  inhomogeneities  [Urick,  1967].    Above 
a  few  hundred  cycles  the  microstructure  becomes  increasingly  important 
and  affects  sonar  transmissions  adversely. 

Sagar  [1960],  by  taking  temperature  profiles  while  collecting  acoustic 
(14.  5  kHz)  transmission  data, was  able  to  associate  acoustic  intensity 
fluctuations  v/ith  the  presence  of  microstructure  while  in  the  absence  of 
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microstructure  the  variation  was  small  and  independent  of  the  Mintzer 
dependence  on  range.  This  illustrates  the  necessity  for  the  Navy  to  be 
concerned  with  thermal  microstructure. 

B.     OPERATIONAL  METHODS. 

XBT  data  is  being  used  by  the  Navy  to  obtain  the  gross  temperature 
profile  for  use  in  the  solution  to  ray  path  equations  and  for  diagrams  to 
predict  sonar  performance.    This  can  be  accomplished  on  board  a  ship 
or  can  be  accomplished  by  Fleet  Numerical  Weather  Central ,  Monterey 
in  a  more  sophisticated  manner,  and  the  products  supplied  as  a  service  to 
fleet  units. 

On  board  ships  the  parameters  most  used  are  the  sea  surface  temper- 
ature and  layer  depth.     Tables  can  then  be  utilized  to  find  the  predicted 
range  for  active  and/or  passive  sonar  contacts.     Little  information  con- 
cerning the  temperature  structure  is  presently  being  used  in  weapon 
selection  or  in  determining  the  expected  performance  of  the  weapon 
once  in  the  environment. 

FNWC ,  Monterey  supplies  information  on  request  to  support  ASW 
operations.    The  products  include  sonar  range  predictions  and  propagation 
loss  predictions.     The  major  factors  which  are  used  for  the  prediction 
computations  are  water  mass  type ,  bottom  type ,  and  depth  of  the  water. 
Temperature  information  is  furnished  from  the  historical  data  bank  with 
no  real  time  input  from  the  area  in  question.    Historical  data  is  used  to 
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prevent  the  computer  solution  from  being  biased  in  case  the  ship  XBT  was 
from  a  small  patch  of  water  rather  than  being  a  representative  sample  of 
the  total  area.    Again,  no  microstructure  information  is  being  used. 

C.     POSSIBLE  UTILIZATION  OF  MICROSTRUCTURE  INFORMATION 
ON  AN  XBT  TRACE . 

A  ship  has  the  capability  of  sampling  the  temperature  in  its  area  when 
desired.    However,  the  XBT  trace  is  not  being  fully  utilized.    By  closer 
examination  of  the  trace ,  the  sonar  operator  could  derive  an  indication 
of  the  quantity  of  microstructure  present  in  his  operating  area.    With  this 
additional  knowledge  the  operator  would  know  whether  to  expect  a  variation 
in  sonar  reception  or  to  expect  steady  reception.     He  would  also  know  of 
the  possibility  of  picking  up  contacts  beyond  the  predicted  range  that 
might  be  held  only  intermittently  if  microstructure  is  present.     He  could 
request  reduced  range  to  maintain  contact  if  the  contact  is  gained  at  the 
predicted  range  and  microstructure  is  present.    Use  could  also  be  made 
of  microstructure  information  in  selecting  the  range  for  launching  a  homing 
torpedo.    The  variation  in  acoustic  intensity  will  be  larger  for  the  f requen- 
cy  range  of  the  homing  torpedo  than  it  will  be  for  sonar  frequencies ,  so 
it  would  enhance  the  opportunity  for  success  of  the  torpedo  if  the  presence 
or  absence  of  thermal  microstructure  is  considered. 

FNWC ,  Monterey  could  aid  the  ASW  problem  by  predicting  the  relative 
magnitude  of  microstructure  to  be  expected  for  an  area  by  examining  the 
microstructure  on  digitized  traces  now  in  the  data  bank  for  the  area.    This 
could  be  made  a  part  of  the  sonar  range  prediction  now  aw  liable. 
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VI.    CONCLUSION 

A.  SUMMARY 

The  motive  behind  the  investigation  of  the  thermal  microstructure 
distribution  in  the  north  Pacific  Ocean  pursued  in  this  thesis  has  been  to 
determine  what  information  can  be  obtained  from  XBT  records  and  then 
used  to  give  some  prediction  of  the  influence  on  sonar  transmissions  and 
receptions.    Various  techniques  were  explored  to  determine  a  measure  of 
the  thermal  microstructure  on  an  XBT  trace  ending  with  the  selection  of 
a  modified  binomial  filter,  suggested  by  Roden  [1971],  being  preferred. 
Data  was  extracted  from  Fleet  Numerical  Weather  Central,  Monterey's 
XBT  data  bank,  for  four  of  Tully's  [1964]  climatic  regions  of  the  north 
Pacific  Ocean. 

B.  CONCLUSIONS 

Based  on  the  examination  of  the  data  from  the  Subtropic ,  Subarctic, 
Kuroshio  Current ,  and  Confluence  regions ,  the  following  conclusions  are 
offered: 

1.  The  modified  binomial  filter  is  a  technique  that  will  produce  a 
mean  temperature  perturbation  vertical  thickness  scale  for  the  thermal 
microstructure  found  on  an  XBT  trace. 

2.  The  digitizing  process  used  by  Fleet  Numerical  Weather  Central, 
Monterey,  divides  the  trace  into  layers  of  linear  temperr  ure  gradient 
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so  that  the  use  of  the  digitized  inflection  points  will  give  the  same  result 
as  that  obtained  using  inflection  points  from  the  temperature  perturbation 
curve. 

3.  The  binomial  filter  will  also  produce  a  measure  of  the  amplitude 
of  the  temperature  perturbations. 

4.  The  average  mean  temperature  perturbation  vertical  thickness 
scale  increases  as  latitude  increases,  while  the  average  maximum  ampli- 
tude of  the  temperature  perturbations  decreases  as  latitude  increases 
in  the  north  Pacific  Ocean  between  20N  and  55N. 

5.  The  average  maximum  amplitude  of  the  temperature  perturbations 
decreases  if  the  average  mean  thickness  scale  is  increased  in  the  north 
Pacific  Oc.pan. 

6.  The  average  mean  temperature  perturbation  vertical  thickness 
scale  was  found  to  be  more  dependent  on  the  latitude  than  Tully's  Climatic 
Regions. 

7.  No  association  could  be  determined  between  the  time  of  day 
or  season  and  the  mean  thickness  scale. 

C.     RECOMMENDATIONS 

The  XBT  digitized  data  could  be  used  for  a  variety  of  studies.    Recom- 
mendations for  further  study  enumerated  below  are  specifically  related  to 
the  development  of  the  data  processed  in  this  thesis. 

1.      One  location  should  be  selected  for  an  extensive  study  to  deter- 
mine if  the  XBT  data  will  reveal  a  correlation  with  time  of  day  and  season. 
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Two  likely  areas  in  the  Pacific  Ocean  for  such  a  study  are  the  San  Diego 
and  Hawaiian  operating  areas  used  by  the  Navy.    XBT  data  should  be  very 
plentiful  from  these  areas. 

2.  A  study  should  be  conducted  to  determine  the  sonar  transmission 
and  reception  variation  that  can  be  expected  from  the  influence  of  thermal 
microstructure  that  appears  on  XBT  traces  taken  during  the  study. 

3.  Fleet  Numerical  Weather  Central,  Monterey,  California,  should 
investigate  the  feasibility  of  including  a  microstructure  quantitative  pre- 
diction as  a  part  of  the  Sonar  Performance  forecast  for  the  active  sonars. 

4.  Fleet  Numerical  Weather  Central,  Monterey,  California,  should 
add  to  the  identification  portion  of  each  digitized  XBT  trace  the  ciimatol- 
ogical  data ,  such  as  wind  speed,  sea  surface  temperature,  and  cloud  cover, 
for  ease  in  evaluation  of  the  trace. 
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